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ABSTRACT 

The C II 133.5nm lines are important observables for the NASA/SMEX mission Interface Region Imaging 
Spectrograph (IRIS). To make 3D non-LTE radiative transfer computationally feasible it is crucial to have a 
model atom with as few levels as possible while retaining the main physical processes. We here develop such 
a model atom and we study the general formation properties of the C II lines. We find that a nine-level model 
atom of C I-C III with the transitions treated assuming complete frequency redistribution (CRD) suffices to 
describe the C II 133.5 nm lines. 3D scattering effects are important for the intensity in the core of the line. The 
lines are formed in the optically thick regime. The core intensity is formed in layers where the temperature is 
about 10 kK at the base of the transition region. The lines are 1.2-4 times wider than the atomic absorption 
profile due to the formation in the optically thick regime. The smaller opacity broadening happens for single 
peak intensity profiles where the chromospheric temperature is low with a steep source function increase into 
the transition region, the larger broadening happens when there is a temperature increase from the photosphere 
to the low chromosphere leading to a local source function maximum and a double peak intensity profile with 
a central reversal. Assuming optically thin formation with the standard coronal approximation leads to several 
errors: Neglecting photoionization severly underestimates the amount of C II at temperatures below 16 kK, 
erroneously shifts the formation from 10 kK to 25 kK and leads to too low intensities. 

Subject headings: line: formation - radiative transfer - Sun: atmosphere - Sun: chromosphere - Sun: transition 
region 


1. INTRODUCTION 

In radiative equilibrium, the temperature naturally de¬ 
creases with height. It therefore came as a surprise when it 
was realised that the solar outer atmosphere shows a temper¬ 
ature increase, from about 6000 K at the visible surface to 
millions of degrees in the corona. It is clear that the energy 
necessary to sustain these high temperatures originates in the 
solar convection zone but the mechanisms of transportation 
and dissipation of the energy to the outer solar atmosphere 
are still hotly debated. 

A central role in this puzzle is played by the solar chromo¬ 
sphere — the region a few thousand kilometres thick between 
the solar photosphere and the transition region and corona. 
This is the region where we go from a magnetic field being 
pushed around by the plasma motions in the lower parts (high 
value of plasma j3, the ratio of the gas pressure to the magnetic 
pressure) to a region where the magnetic pressure dominates. 
The energy that heats the solar corona needs to be transported 
through the chromosphere. The total energy needed to bal¬ 
ance the radiative losses in the chromosphere is also at least 
ten times larger than what is needed to sustain a hot corona 
(Withbroe & Noyes 1977). 

In order to spectroscopically diagnose the chromosphere, 
we need spectral lines with high enough opacity to have the 
formation there. This means resonance lines of abundant el¬ 
ements from the ionization state that dominates under chro¬ 
mospheric conditions or lines from excited levels of hydro¬ 
gen or helium. In the optical part of the spectrum accessible 
from ground based observatories we are restricted to a hand¬ 
ful of spectral lines that meet these criteria: the most impor- 
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tant being the hydrogen H-a line, the resonance lines and in¬ 
frared triplet lines from singly ionized calcium and the helium 
1083 nm line. These diagnostic lines have provided a wealth 
of information on the solar chromosphere but there are chal¬ 
lenges in the interpretation of the observations: the H-a line 
has a very large thermal width and is a very strongly scattering 
line largely decoupled from the local conditions (Leenaarts 
et al. 2012), the calcium lines carry the chromospheric sig¬ 
nal only in the core of a strong photospheric absorption line 
and the interpretation is heavily dependent on numerical mod¬ 
elling (Carlsson & Stein 1997) while the helium 1083 nm line 
has a complicated formation with a lower level mainly popu¬ 
lated from ionization from coronal radiation followed by re¬ 
combination. 

The ultraviolet (UV) part of the spectrum provides many 
more diagnostic possibilities for chromospheric studies. 
There are many resonance lines from abundant elements in 
this part of the spectrum, the continuum radiation is much 
fainter in the UV, and the chromospheric diagnostics are in¬ 
trinsically stronger. The recently launched NASA/SMEX 
mission Interface Region Imaging Spectrograph (IRIS) was 
especially designed to target the solar chromosphere and tran¬ 
sition region (De Pontieu et al. 2014). IRIS provides high 
spectral, spatial and temporal resolution in three wavelength 
bands: 133.2— 135.8nm and 138.9 — 140.7nm with 1.3 pm 
wavelength pixels (corresponding to 2.9 — 2.7 km s -1 ) and 
278.3 — 283.4nm with 2.5 pm wavelength pixels (2.7 km s -1 ), 
0.167 arcsec spatial pixels and 0.33 arcsec slit width. 

The resonance lines from singly ionized carbon (C II) 
around 133.5 nm are among the strongest lines in the solar 
ultraviolet spectrum. There are three components in the multi- 
plet: at 133.4532nm and at 133.5708nm with a weaker blend 
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at 133.5663 nm. The large number of photons emitted in the 
lines makes it possible to observe at high spatial and tempo¬ 
ral resolution and the C II lines are therefore among the more 
important diagnostic lines for IRIS. To explore this diagnostic 
potential it is important to find out how the lines are formed 
and thereby establish how the conditions in the solar outer 
atmosphere are encoded in the line profiles. Singly ionized 
carbon is the dominant ionization state throughout most of the 
upper chromosphere but C II also exists in significant amounts 
up to middle transition region temperatures, as we will show 
here. 

The very first C II solar spectrum was taken by the normal- 
incidence grating spectrograph from the U.S. Navel Research 
Laboratory (NRL) April 19, 1960 (Detwiler et al. 1961). 
Later, the lines have been observed several times from the 
series of Orbiting Solar Observatory (OSO) satellites, espe¬ 
cially OSO-4, although this instrument only obtained barely 
resolved C II spectra (Pottasch 1964). These authors sug¬ 
gested that the primary emission comes from regions with 
a temperature of 20 kK. Later, Chipman (1971) used OSO- 
4 and OSO-6 observations of the C II lines to constrain the 
temperature in the upper chromosphere and found the obser¬ 
vations to be compatible with a temperature plateau at about 
15 kK. Lites et al. (1978) used OSO-8 observations to con¬ 
strain the temperature structure of the upper chromosphere 
using centre-to-limb measurements of the C II lines. Through 
spectrum synthesis of the lines they found them to be ex¬ 
tremely sensitive to the temperature and physical extent of the 
plateau at 20 kK proposed by Vernazza et al. (1973) and found 
the observations more compatible with a plateau at 16.5 kK 
with 25% more material than in the models by Vernazza et al. 
(1973). The OSO-8 observations reported on by Lites et al. 
(1978) were the highest spectral resolution observations thus 
far and the average quiet Sun profiles show a clear central re¬ 
versal and a red emission peak slightly brighter than the blue 
peak. 

A detailed discussion of the formation of the C II 133.5nm 
multiplet appears in Avrett & Loeser (2008). This study 
is based on a one-dimensional, time-independent semi- 
empirical model constructed to reproduce the average quiet- 
Sun observations presented in the SUMER atlas of the ex¬ 
treme ultraviolet spectrum (Curdt et al. 2001). Avrett et al. 
(2013) presented further calculations of the C II 133.5nm 
multiplet in four semiempirical models representing the faint 
and mean internetwork, a network lane, and bright network. 

The C II lines have also been used to diagnose other stars. 
Brown & Carpenter (1984) used them to constrain the chro¬ 
mospheric temperature structure in late type giant and super¬ 
giant stars. The C II multiplet is also considered as a good 
diagnostic to estimate stellar basal fluxes (Schrijver 1995). 
A comprehensive study based on the analysis of time series 
data with the Solar Ultraviolet Measurement of Emitted Ra¬ 
diation (SUMER) instrument on-board SOHO was presented 
by Judge et al. (2003). By comparing the measured intensi¬ 
ties with synthetic data from simulations of acoustic waves 
they conclude that even the lowest observed levels need to 
have a magnetic heating component in addition to heating by 
the dissipation of acoustic waves, thus questioning the com¬ 
mon thought that the basal flux from stellar chromospheres is 
accounted for by non-magnetic heating. 

With the advent of detailed, 3D radiation magnetohydrody¬ 
namic (RMHD) models spanning from the convection zone 
to the corona it has become possible to study the formation of 
spectral diagnostics in an inhomogeneous and dynamic set¬ 


ting. Such studies provide a more extensive understanding 
of how the atmospheric properties are encoded in the detailed 
line-profiles than the earlier studies based on one-dimensional 
semi-empirical models. In the first papers in this series of pa¬ 
pers on the formation of IRIS diagnostics we have employed 
this approach to study the Mg II h & k lines (Leenaarts et al. 
2013a,b; Pereira et al. 2013) and the Mg II triplet lines (Pereira 
et al. 2015). We will now present results for the C II lines. To 
make 3D non-LTE radiative transfer computationally feasible, 
it is crucial to have a model atom with as few levels as pos¬ 
sible while retaining the main physical processes. In this first 
paper we develop such a model atom and we study the gen¬ 
eral formation properties of the intensity profiles of the C II 
lines. In the next paper we will present statistical correlations 
between the atmospheric parameters and observables. In the 
third paper on the C II lines we will use IRIS observations to 
further test the diagnostic potential of the C II lines and use the 
comparison of observations with the synthetic observables to 
draw conclusions on what might be missing physical ingredi¬ 
ents in the simulations. 

The layout of this paper is as follows: in Section 2 we give 
details on how we solve the coupled equations of statistical 
equilibrium and radiative transfer, in Section 3 we describe 
the various model atmospheres we use, in Section 4 we de¬ 
scribe the model atom, the atomic data and the procedure that 
ends up with a quintessential model atom. In Section 5 we 
present the general formation properties of the C II lines, in 
Section 6 we discuss the formation in a Sunspot atmosphere 
and we end with conclusions and discussion in Section 7. 

2. RADIATIVE TRANSFER 

For the detailed analysis of the C II lines we solve the sta¬ 
tistical equilibrium equations in non-LTE using three codes. 
For the study of our atomic model, basic radiative trans¬ 
fer processes and the simplifications of the model atom we 
use the ID code MULTI (Carlsson 1986). In order to 
study 3D effects we use the full 3D radiative transfer code 
MULTI3D (Leenaarts & Carlsson 2009). To study the im¬ 
portance of partial redistribution (PRD) effects and effects 
of calculating background continua with several elements in 
non-LTE simultaneously we use the code RH (Uitenbroek 
2001). MULTI solves the non-LTE radiative transfer problem 
in semi-infinite, plane parallel one dimensional atmospheres 
with prescribed macroscopic velocity fields. Complete fre¬ 
quency redistribution (CRD) is assumed for all bound-bound 
transitions and the background scattering is considered to be 
coherent. Line transitions that overlap in frequency are not 
treated self consistently (which is of relevance here since one 
of the C II lines is a blend with two components). Version 2.3 
of MULTI includes the local approximate operator of (Olson 
et al. 1986), Ng acceleration (Ng 1974), collisional-radiative 
switching (Hummer & Voels 1988) and the opacities are taken 
from the Uppsala Opacity Package (Gustafsson 1973). The 
background continuum opacity is normally calculated assum¬ 
ing LTE for the relevant atoms and ions but it is possible to 
iteratively include non-LTE populations for the background 
opacity. 

The full three dimensional radiative transfer code 
MULTI3D is not only superior in that it treats the non-LTE 
problem in full 3D but the treatment of overlapping lines and 
continua is self-consistent and not based on iteration. How¬ 
ever, the computations are naturally more computationally in¬ 
tensive in 3D which is why we use the ID version for experi¬ 
mentation with large model atoms. We assume CRD also for 
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the computations with MULTI3D. 

RH is another non-LTE code capable of treating PRD in the 
lines and calculating overlapping continua by solving for sev¬ 
eral elements in non-LTE simultaneously. This code is there¬ 
fore used for the tests of the significance of PRD and for the 
detailed continuum calculations. 

3. MODEL ATMOSPHERES 

We use a number of different model atmospheres in this 
work. We start by studying the atomic processes important 
for the lines we are interested in by using the often used 
semi-empirical VAL3C model atmosphere (Vernazza et al. 
1981). Once we have arrived at a simplified atomic model that 
preserves the important characteristics of the more extended 
atomic model, we employ a 3D model atmosphere to test the 
simplified model atom over a wider parameter range and to 
see the sensitivity of the line formation to atmospheric prop¬ 
erties. To study radiative transfer effects in 3D, we employ the 
full 3D solution. The 3D atmosphere is taken from a 3D sim¬ 
ulation with the Bifrost code (Gudiksen et al. 2011). Bifrost 
solves the equations of radiation magnetohydrodynamics on 
a staggered Cartesian grid. Detailed radiative transfer is 
included through the multi-group opacity approach (Nord- 
lund 1982) modified to take into account scattering (Skartlien 
2000). Radiative losses in the chromosphere are calculated 
in non-LTE using simplified recipes (Carlsson & Leenaarts 
2012) that are based on detailed ID full non-LTE radiative 
transfer simulations. Non-equilibrium ionization of hydro¬ 
gen is included following the description of Leenaarts et al. 
(2007). The simulation we use encompasses 24 x 24 x 16.8 
Mm on the Sun discretised onto a 504 x 504 x 496 grid. Ver¬ 
tically, the computational volume extends from 2.4 Mm be¬ 
low to 14.4 Mm above average optical depth unity at 500 nm 
(which is the zero-point of our height-scale), and thus covers 
the solar atmosphere from the upper convection zone, photo¬ 
sphere, chromosphere to the lower corona. Horizontally, the 
grid is equidistant with a spacing of 48 km. The z-axis is not 
equidistant. It has grid spacing of 19 km between z = — 1 and 
z = 5 Mm, whereas the spacing increases towards the lower 
and upper boundaries to a maximum of 98 km at the coro¬ 
nal boundary. Magnetic fields are introduced into a relaxed 
hydrodynamic simulation by specifying the vertical magnetic 
field at the bottom boundary with a potential field extrapo¬ 
lation producing the magnetic field throughout the 3D box. 
In our particular simulation we specify two concentrations of 
opposite polarities separated by 8 Mm. The potential field is 
quickly swept to the intergranular lanes and slowly builds up 
a set of loops between the opposite polarities. We have cho¬ 
sen a snapshot at a time 3850s after the magnetic field was 
introduced, well after initial transients have passed through 
the box and the heating of the atmosphere through the work 
of the convection on the magnetic field has come to a quasi¬ 
steady state. This simulation is described in detail in Carlsson 
et al. (2015) and is the same as used in Leenaarts et al. (2012), 
Stepan et al. (2012), de la Cruz Rodriguez et al. (2013) and 
for the first papers in this series of papers on the formation of 
IRIS diagnostics (Leenaarts et al. 2013a,b; Pereira et al. 2013, 
2015). 

The vertical magnetic field strength in the photosphere is 
shown in Figure 1 . The two opposite magnetic field polarity 
patches are clearly seen. We will illustrate formation prop¬ 
erties along a slice through the model at v= 12 Mm, marked 
with a vertical line, and at four different positions, marked 
with letters A-D in Figure 1 . 


To study the formation of the C II lines also in a more ac¬ 
tive atmosphere, we employ a semi-empirical atmosphere of 
a sunspot umbra. The standard model SPOTM of Maltby 
et al. (1986) has an emphasis on the temperature structure of 
the photosphere up to the temperature minimum region and 
does not extend high enough to cover the formation region 
of the C II lines. We therefore use the semi-empirical model 
of Lites & Skumanich (1982) that is based on OSO-8 obser¬ 
vations of the chromospheric lines of hydrogen (Ly-a, Ly-/3), 
Mg II (k,h), Ca II (K,H) and the transition region C IV 1548 nm 
resonance line. 
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Figure 1. Vertical magnetic field of the Bifrost atmosphere at z = 0 height. 
The vertical line shows the position of the 2D cut used and the letters A-D 
mark columns chosen for detailed study in section 5.1. 

The temperature structures of the two semi-empirical mod¬ 
els and the four columns in the Bifrost snapshot are shown 
in Figure 2. Both the semi-empirical models have an artificial 
temperature plateau at 24 kK in order to get enough flux in the 
hydrogen Lyman-a line. The umbral model has this plateau 
and the transition region at a slightly smaller column mass 
than the quiet Sun model VAL3C. The transition region in the 
Bifrost model is typically at an even smaller column mass. 

4. QUINTESSENTIAL MODEL ATOM 

One of the objectives of the current work is to arrive at a 
simplified model atom that includes the non-LTE processes 
important for the formation of the intensity profiles of the 
C II 133.5 nm lines but is small enough to enable large scale 
3D non-LTE simulations. To that end we start with a more 
complex model atom, the one used by Judge et al. (2003). 
This atomic model contains a total of 22 levels: 8 levels of 
C I, 6 levels of C II, 4 levels of C III, 3 levels of C IV and 
the ground state of C V. We simplify this model by remov¬ 
ing and merging levels following the procedure described in 
Bard & Carlsson (2008). At each step we solve the non-LTE 
problem and compare the intensities in the C II 133.5 nm lines 
and the ionization fraction as function of depth with that of 
the full model atom. As a first step we remove all levels 
above the ground state of C III. This results in a model atom 
of 15 levels. In the next step we remove levels not impor¬ 
tant for the ionization balance and the levels not directly cou¬ 
pled to the lower and upper levels of the C II 133.5 nm lines. 
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Figure 2. Temperature as function of logarithmic column mass for atmo¬ 
spheric models in this paper: The quiet Sun model VAL3C (thick grey), the 
umbral model of Lites & Skumanich (1982) (thick black), and four columns 
from the Bifrost snapshot as marked in Figure 1: A (blue), B (red), C (brown) 
and D (green). 

The quintessential model atom we arrive at has a total of 9 
levels, see Figure 3. The lines we are interested in are the 
unblended component at 133.4532nm (hereafter called the 
133.4nm line) and the stronger 133.5708nm line that has a 
weaker blend at 133.5663nm, hereafter collectively referred 
to as the 133.5nm line. A comparison between the results 
for the VAL3C atmosphere based on the 22-level model atom 
and the 15 and 9-level atoms is given in Figure 4. Both the 
intensities and ionization structure are almost identical calcu¬ 
lated from the quintessential 9-level model atom as from the 
full 22-level atom. This is also true for the full Bifrost atmo¬ 
sphere — on average the quintessential 9-level model has 5% 
larger peak intensity in the C II lines but less than 5% of the 
columns show a difference of more than 15% and the largest 
difference is 27% compared with the full 22-level atom. 

We have also made experiments with a more complete C I 
model atom including highly excited levels to get recombina¬ 
tion channels right and including the photoionizing radiation 
from the hydrogen Lyman-a line. This influences the C I/C II 
balance at the low temperature end but has negligible influ¬ 
ence on the CII 133.5 nm lines, see Section 4.3. 

We stress that our model atom has been constructed for the 
specific purpose of modeling the intensity of the CII 133.5 nm 
lines. If other lines are of interest or if polarization is to be cal¬ 
culated, this model atom might need to be extended to include 
additional levels. 


Figure 3. C II atomic model with 9 levels. The C II multiplet of special 
interest here consists of the three allowed transitions between the 2s2p 22 D 
and the 2 s 1 2 2p 2 P° terms. 



Figure 4. Intensities and ionization fractions resulting from the original 
22-level atom and the simplified 15-level atom (panels (a) and (c)) and the 
quintessential 9-level atom (panels (b) and (d)). Results for the reference 
model are shown in black and for the 15/9-level models in red. The VAL3C 
model atmosphere was used and the weak blend at 133.566 nm was neglected. 


4.1. Atomic data 

In the 9-level atomic model level energies for C I are from 
the NIST 1 online database, for C II from Moore (1993) and the 
C III level is from Martin et al. (1978). The transition proba¬ 
bilities are from the NIST database. Detailed photoionization 
data (from C I to C II) are from the TOPBASE 2 database of 
the OPACITY project. The bound-bound collisional rates for 
C II are from Tayal (2008). The bound-bound collisional rates 
for neutral carbon are determined using the impact parameter 
method of Seaton (1962). As we will see later, the intensity 
ratio between the multiplet components depends on the rela¬ 
tive population of the upper levels. As the energy difference 

1 http://aeldata.nist.gov/nist_atomic_spectra.html 

2 http://cdsarc.u-strasbg.fr 


between the levels of the 2 p 2 D term is only 2.5 cm -1 (0.008 
eV), collisions with neutral hydrogen may be important. For 
these collisions we use data from Bahcall & Wolf (1968) with 
P(i -A f) = 4/10 (Barklem 2014, private communication). 

In our simulations the total recombination rates include 
the radiative rates, dielectronic recombination rates and three 
body recombination rates. The radiative recombination rates 
are either calculated from the detailed photoionization cross- 
sections (recombination into C I) or as tabulated by Shull & 
van Steenberg (1982). The three body recombination rates 
are from Arnaud & Rothenflug (1985) and the dielectronic 
recombination rates are from Altun et al. (2004) (C II -A C I) 
and Colgan et al. (2003) (C III -A C II). We set the carbon 
abundance to 8.43 (Asplund et al. 2009) on the usual logarith¬ 
mic scale where the abundance of hydrogen is 12. 
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4.2. Continuum intensity 

In addition to the opacity from the multiplet lines we are in¬ 
terested in, we have a background opacity that varies slowly 
with wavelength. At 133.5 nm under solar chromospheric 
conditions, this background opacity is dominated by bound- 
free transitions in silicon, carbon, hydrogen and iron with 
electron scattering taking over as the dominant background 
opacity source in the upper chromosphere and corona, see 
Figure 5. The carbon background opacity comes from bound- 
free transitions from the levels 2 p 2 l D and 2 p 2 l S while the 
ground term 2 p 2 3 P has a bound-free edge at 1 lOnm and does 
not contribute. Note that the energy difference from the S I 
levels 3 p 4 l D and 3 p 4 l S and the ground state of SII places the 
bound-free edges longward of our carbon multiplet and these 
transitions would then contribute to the background opacity. 
However, these levels ionize to an excited state of S II placing 
the bound-free edges shortward of our multiplet and sulphur 
therefore does not contribute significantly to the background 
opacity. Hydrogen, carbon and silicon are all out of LTE in 
the chromosphere so a proper calculation of the continuum 
intensity needs to treat all three elements in non-LTE simul¬ 
taneously. Our carbon model atom includes the singlet levels 
of C I so that contribution is in non-LTE automatically. Hy¬ 
drogen non-LTE level populations are also provided in the at¬ 
mospheric models used here so the only remaining problem 
is the silicon contribution. To include silicon in non-LTE si¬ 
multaneously with carbon would be computationally costly. 
This is of course necessary if one is interested in getting the 
continuum intensity correct. However, in cases where the con¬ 
tinuum intensity is very much lower than the carbon multiplet 
intensity, the correct treatment of the background continuum 
may be unimportant for getting the correct multiplet intensity. 
We tested this by comparing the full treatment with one where 
silicon is assumed to be in LTE for the calculation of the back¬ 
ground opacity, see Figure 6. It is clear that the continuum 
intensity is severely overestimated when silicon is treated in 
LTE but that such a simplified treatment does not significantly 
affect the line intensity. 



Figure 5. Relative contribution to the background opacity at 133.5nm as a 
function of height in the VAL3C atmosphere. The source of the opacity is 
bound-free opacity from Si, C, Fe and H and electron scattering ( o(e ~)) and 
Rayleigh scattering on neutral hydrogen (R(H)). The secondary height scale 
shown in red is the logarithm of the continuum optical depth at 133.5nm. 
Hydrogen, carbon and silicon have all been treated in non-LTE while iron 
was treated in LTE. 




133.43 133.44 133.45 133.46 133.47 133.54 133.55 133.56 133.57 133.58 133.59 

Wavelength [nm] Wavelength [nm] 

Figure 6. C II line intensity for the C II 133.4nm (left) and C II 133.5nm 
(right) lines, when assuming non-LTE Si background opacity (red) and LTE 
Si background opacity (black). The VAL3C atmosphere was used for the 
comparison. Note the logarithmic intensity scale. 

4.3. Ionization balance 

The formation of the C II-lines is to a large extent set by 
the ionization balance - where in the atmosphere is the C II 
ionization stage dominant? This in turn is set by the ioniza¬ 
tion and recombination rates between C II and C III and by 
the rates between C I and C II. The balance between C I and 
C II is dominated by photoionization and radiative recombi¬ 
nation. To get this balance right, we must include the levels 
that dominate the photoionization (normally the ground lev¬ 
els and levels with low excitation) and highly excited levels 
that dominate in the recombination channels. The photoion¬ 
ization of C I is dominated by ionization from the 2 p 2 3 P term 
with edges around 110 nm and from the 2 p 2l D level with 
photoionization edges close to 123.9 nm. For the latter, pho¬ 
toionization from hydrogen Lyman-a radiation may play an 
important role. This complicates the computational task enor¬ 
mously since the proper solution involves a large number of 
C I energy levels and the simultaneous solution of the hydro¬ 
gen non-LTE problem with the Lyman-a line treated in partial 
frequency redistribution. We have made test calculations with 
such a setup and find that by neglecting the photoionization 
from Lyman-a and all levels above the ground term and the 
two first excited levels of C I, the amount of C II at tempera¬ 
tures below 10 kK is indeed underestimated but the effect on 
the C II lines is minimal (intensity changes by less than 5%). 

In the classical coronal ionization equilibrium approxima¬ 
tion, photoionization is neglected and the ionization rate is 
assumed to be dominated by collisional ionization. This leads 
to a severe underestimate of the amount of C II below a tem¬ 
perature of 20 kK under solar chromospheric conditions, as is 
clearly shown in Figure 7. 

At the high temperature end, the ionization fraction of C II 
is determined by collisional ionization and dielectronic re¬ 
combination. Part of the dielectronic recombination is to high 
energy levels and there may be collisional ionization taking 
place before the electron cascades down to lower energy lev¬ 
els through allowed radiative transitions. A correction in the 
dielectronic recombination rates is made following Summers 
(1972) as implemented by Judge (2007) by including a den¬ 
sity dependent factor in the dielectronic recombination coef¬ 
ficient. This decreases the recombination rate and we get a 
smaller C II-fraction at the high temperature range than if the 
dielectronic recombinations to high energy C II-levels are all 
assumed to end up with radiative cascades to the ground state, 
see Figure 7. Neglecting Summer’s correction factor for di¬ 
electronic recombination, the C II ionization fraction does not 
fall below 10% until a temperature of 75 kK is reached. In¬ 
cluding this factor, the fraction of C II drops rapidly at a tern- 
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perature of 25 kK and we have an ionization fraction below 
10% already at 50 kK. 



T [kK] 

Figure 7. CII ionization fraction in the VAL3C model atmosphere using the 
adopted atomic parameters (solid), neglecting photoionization from CI to CII 
and also neglecting density dependent dielectronic recombination (=default 
Chianti ionization balance) (dashed). Note the logarithmic temperature scale. 

The radiative and dielectronic recombination rates are a 
function of temperature and linearly dependent on electron 
density apart from the above mentioned density correction of 
the dielectronic recombination rate. Therefore, if the ioniza¬ 
tion is dominated by collisional ionization (also a function of 
temperature and linear in electron density) and three body re¬ 
combination is negligible, the ionization balance is a function 
of temperature only. In Figure 8 we show the fraction of C II 
in the full Bifrost snapshot (504x504x496 grid points). It is 
clear that in the high temperature end we are close to having 
an ionization balance only dependent on temperature. The 
spread is caused by three-body recombination and the density 
dependent dielectric recombination. At the low temperature 
end, the ionization is dominated by photoionization and the 
ionization fraction is not a single valued function of tempera¬ 
ture. There is also an extra spread caused by the 3D radiation 
field. 



Figure 8. Probability density function of the CII ionization fraction as func¬ 
tion of temperature in the Bi frost simulation. In each temperature bin the 
PDF has been normalized to the largest value for increased visibility. Note 
the logarithmic temperature scale. 


4.4. Effect of Partial frequency Redistribution 

The C II lines are resonance lines (lower levels part of 
the ground term) and one might expect effects of partial fre¬ 
quency redistribution of photons (PRD) to be important. We 
tested this by solving the non-LTE problem for given atmo¬ 
spheres both including PRD effects and using the computa¬ 
tionally advantageous approximation of complete redistribu¬ 
tion (CRD). We did not find any significant differences, see 
Figure 9, and conclude that CRD can be used in the modelling 
of the C II 133.5 nm lines under solar conditions. 



Figure 9. Comparison between partial redistribution (PRD, red) and com¬ 
plete redistribution (CRD, black) approximation for the C II lines in the 
VAL3C atmosphere for the C II i33.4nm line (a) and C II 133.5nm blend 
(b). 


4.5. Importance of 3D radiative transfer 

Figure 10 shows a comparison of the vertically emergent 
intensity in the core of the C II 133.4nm line solving each 
column in the 3D atmosphere independently as a ID problem 
and solving the full 3D radiative transfer problem. It is clear 
that 3D scattering effects are important in the core of the line 
and a full 3D transfer solution should be employed. It also 
means that we should expect the intensity close to the core 
to be influenced by the surroundings and not provide a very 
accurate mapping of the conditions along the local column. 



0 5 10 15 20 0 5 10 15 20 

x [Mm] x [Mm] 


Figure 10. Comparison of the vertically emergent intensity in the core of the 
C II 133.4nm line solving each column in the 3D atmosphere independently 
as a ID problem (left) and solving the full 3D radiative transfer problem 
(right). 

5. RESULTS 

In this section we will study in detail how the C II lines at 
133.5 nm are formed. We use the quintessential 9-level model 
atom described above, assuming CRD but using the full 3D 
solution for the Bifrost atmosphere results. 
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5.1. Contribution function to the intensity 

As a tool to find out which parts of the atmosphere con¬ 
tribute to the emergent intensity of the C II lines, we use the 
contribution function to intensity. 

The emergent intensity in a ID plane parallel semi-infinite 
atmosphere can be written as 

Iv,n= T-Sy(z)e~^^Xvdz ( 1 ) 

io M 

In this equation, /I = cos0, where 6 is the angle between 
the line of sight to the observer and the local vertical. The 
source function (5 V ), opacity (% y ) and optical depth (t v ) are 
functions of frequency (v) and geometrical height (z). At disk 
center (jU = 1), the integrand in Equation 1 describes the local 
creation of photons (, SyXvdz ) and the fraction of those that es¬ 
cape (e~ Zy ). It is thus natural to define the contribution func¬ 
tion to intensity on a geometrical height scale as 

Ci v (z) = S v e-*Xv (2) 

Following Carlsson & Stein (1997) we also rewrite the con¬ 
tribution function as 

C Iv (z)=S v z v e- Zv ^. (3) 

T v 

where the term z v e~ Zy having a maximum at T v = 1 repre¬ 
sents the Eddington-Barbier part of the contribution function, 
S v gives the source function contribution and the final term, 
^ picks out effects of velocity gradients in the atmosphere. 

To illustrate the range of line formation scenarios we have 
picked four different columns from the Bifrost atmosphere. 
These columns are marked in Figure 1 . 

Figures 11-14 show the total contribution function as func¬ 
tion of height and frequency across the CII 133.4nm line (bot¬ 
tom right panels) as well as the three individual terms above 
in the same manner as in the four-panel diagrams in Carlsson 
& Stein (1997) for these four columns. To avoid the confusion 
from the blend we have chosen the 133.4nm line to analyze 
the formation in detail but the same general characteristics 
are true for the stronger line as well but with a formation a bit 
higher in the atmosphere. 

Figure 11 shows a typical case of a double peaked profile 
(the intensity profile is shown in the bottom right panel). The 
velocity (turquoise almost vertical line in all panels) does not 
exhibit strong gradients which is why the % v / T v term in the 
upper left panel does not exhibit strong variations with wave¬ 
length or height. The T v = 1 (red dashed line in all panels) 
ranges from 0.8 Mm in the continuum to 2.1 Mm at line cen¬ 
ter. The source function (solid green line in the upper right 
panel) is equal to the Planck function (the dashed yellow 
line shows the radiation temperature of the Planck function 
and thus the temperature itself) from the photosphere up to 
a height of 0.7 Mm where it decouples. The strong peak in 
temperature at a height of 1.6 Mm causes a peak in the source 
function (not so pronounced in the radiation temperature rep¬ 
resentation shown with the green line but quite evident in the 
linear representation shown in the image in the upper right 
panel). This peak in the source function is responsible for the 
two emission peaks in the intensity profile. Because of the 
lack of strong velocity gradients, the intensity profile is quite 
symmetric with almost equal intensity of the read and blue 
peaks. The line core is formed at the highest height of the op¬ 
tical depth unity curve at 2.1 Mm. The source function is there 


lower than at 1.6 Mm and we get a central reversal. At 2.1 Mm 
height there is a downward velocity of 2 km s -1 which causes 
a redshift of the line center with a similar amount. 

Figure 12 shows a case of a single peaked profile where the 
formation is optically thin. It is clear from the total contribu¬ 
tion function image (lower right panel) that the maximum of 
the contribution to intensity comes from a height of around 
2.6 Mm which is well above the maximum height of optical 
depth unity (2.3 Mm). The atmosphere is moving upwards 
with 7 km s -1 where the intensity is formed and this gives a 
blueshift of the same amount. Note that the whole profile has 
a blueshift of 5-7 km s -1 even though the atmospheric veloc¬ 
ity is almost zero below a height of 1.3 Mm. This is because 
most of the intensity profile is formed around 2.6 Mm height 
and not at the optical depth unity height. The asymmetry with 
a stronger red wing than blue wing comes from the velocity 
gradient giving a larger Xv/^v term (upper left panel) on the 
red side. The full-width-half-maximum (FWHM) of the in¬ 
tensity profile is 10.6 km s -1 corresponding to the thermal 
broadening at 29 kK which is consistent with the temperature 
in the lineforming region. 

Figure 13 shows the formation of a single peaked profile 
in a flux concentration (column C in Figure 1 at (x,y) = 

(6.9.12.9) Mm). Because of the strong magnetic field (-1.6 
kG in the photosphere), the density is lower than in the sur¬ 
roundings and a given optical depth occurs at a lower geomet¬ 
rical height. The continuum is therefore formed at 0.25 Mm 
compared with the average continuum formation height of 
0.8 Mm. The core of the line is also formed slightly lower 
than in the surroundings, at 1.9 Mm height, but close to the 
height where the monochromatic optical depth is unity. The 
line thus shows an optically thick formation with the inten¬ 
sity being a map of the source function at optical depth unity. 
The source function is monotonically increasing with height 
up to where the core is formed so there is no central depres¬ 
sion. Note also the rather box-shaped profile - a consequence 
of a source function that increases more rapidly just above the 
continuum formation height than further up towards the core 
formation height. 

Finally, Figure 14 shows the formation of the intensity in 
another flux concentration (column D in Figure 1 at (x,y) = 

(14.5.9.9) Mm). The continuum is formed low because of 
the strong magnetic field (1.7 kG in the photosphere). The 
continuum intensity is high because of the high temperature 
at optical depth unity in the continuum (see also Figure 2). 
Above the photosphere we have strong velocity gradients that 
cause a broadened asymmetric profile. The source function is 
increasing up to a height of 0.6 Mm, is decreasing up to the 
height of optical depth unity at line center and then shows a 
sharp increase.We therefore get a profile with three peaks with 
the outer ones formed at the wavelengths where optical depth 
unity is at 0.6 Mm. The FWHM of the intensity profile is 26.1 
km s -1 which is much larger than the thermal width in the line 
forming region. This is partly because of the large velocity 
gradients in the atmosphere at this column but also the profile 
calculated without velocity gradients in the atmosphere has a 
profile much broader than what would be expected from the 
thermal width. This is typical for a line that is formed in the 
optically thick regime. The line-width depends strongly on 
how the source function varies with height in the atmosphere. 
We will return to this point in Section 5.2 

So far we have exemplified the formation of the 
C II 133.4nm line through a detailed analysis of four differ- 
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Figure 11. Contribution function to intensity of the CII 133.4nm line. Each 
panel shows the quantity specified in the top-left corner as a grey-scale image 
as function of frequency from line center (in Doppler shift units) and atmo¬ 
spheric height z, see equation 3. Multiplication of the first three produces the 
contribution function to intensity shown in the lower right panel. A t v = 1 
curve (red dashed) and the vertical velocity (turquoise solid, positive is up- 
flow) are shown in each panel, with a v z = 0 line in the upper left panel for 
reference. The upper-right panel also contains the Planck function (yellow 
dashed) and the total source function (green solid) in radiation temperature 
units specified along the top. The lower-right panel also contains the emer¬ 
gent intensity profile (white solid), with the scale along the right-hand side. 
The formation diagram is shown for column (jc,y) = (12,10.3) Mm in the 
Bifrost atmosphere (column (A) in Figure 1). 


T[kK] 
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Figure 12. As Figure 11 but for column (x,y) = (12,22.8) Mm in the Bifrost 
atmosphere (column (B) in Figure 1). 


Figure 13. As Figure 11 but for column (x,y) = (6.9,12.9) Mm in the 
Bi frost atmosphere (column (C) in Figure 1). 
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Figure 14. As Figure 11 but for column (x,y) = (14.5,9.9) Mm in the 
Bifrost atmosphere (column (D) in Figure 1). 


ent columns in the Bi frost atmosphere. We will now look at 
the formation of the intensity along the cut at x= 12 Mm (see 
Figure 1). We will investigate the contribution function to 
intensity at the line core; where for each column in the simu¬ 
lation atmosphere, the line core is defined as the wavelength 
with maximum intensity (single peaked profiles) or minimum 
intensity in the central reversal (double peaked profiles). This 
contribution function is given in Equation 2. We are also in¬ 
terested in where the integrated intensity is coming from. We 
achieve this by defining the contribution function to total in¬ 
tensity as the integral over frequency of the contribution func¬ 
tion to specific intensity: 
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Ci t 0 Jz)~Jc Iv (z)dv (4) 

Note that this definition only makes sense when the contin¬ 
uum intensity is much lower than the line intensity (which 
is the case here). We perform the integration over a wide 
enough frequency range to cover the width of the line but nar¬ 
row enough not to be dominated by the continuum intensity 
contribution. 

The top panel of Figure 15 shows the contribution func¬ 
tion to the intensity at the line core along this cut. The 
line core is normally formed between 10 kK and 40 kK but 
there are also regions where the line core contribution func¬ 
tion reaches its maximum at lower temperatures (e.g. around 
y = 17 — 18 Mm). This happens when the transition region 
is at a low density. Optical depth unity is reached when we 
have high enough density. With a low density in the transi¬ 
tion region, this happens deeper down where the temperature 
is lower. Because of the lower temperature, the line is also 
weaker (since the intensity is the integral of the contribution 
function, Figure 15 also gives an indication of the intensity). 
The formation layer of the line core is between 2 and 4 Mm 
height and is as corrugated as the temperature structure. 

The bottom panel of Figure 15 shows the contribution func¬ 
tion to total intensity. The total intensity is mostly coming 
from the core formation region but there is also a broad tail of 
the contribution function going all the way down to the for¬ 
mation region of the continuum (which is on average around 
a height of 0.8 Mm). At y = 0 — 0.5 Mm the contribution func¬ 
tion to total intensity is dominated by the continuum contri¬ 
bution. Here the continuum is enhanced because of a high 
temperature(up to 7 kK) around 0.8 Mm and the line is actu¬ 
ally in absorption. 

The height of unit optical depth (and therefore also the for¬ 
mation height) varies with the density stratification. Since 
the 133.5nm lines are resonance lines, we expect the optical 
depth unity to be at given column density of singly ionized 
carbon. The opacity at the line core depends on the broad¬ 
ening of the profile but is on the order of 10 6 cm 2 g _1 for the 
C II 133.4 nm line. We therefore expect optical depth unity 
to be at a column mass of about 10 -6 gcm -2 . We can define 
a contribution function on a logarithmic column mass scale 
from 

C/ V (log 10 m c )d(log 10 m c ) = Q v (z)dz (5) 

where m c is the column mass, which leads to 

O v (l°gio w c) = ^ln(10 )C/ V (z). (6) 

Figure 16 shows the contribution function to intensity at 
line core (top panel) as well as the contribution function to 
total intensity (bottom panel) on a logarithmic column mass 
scale. The range in logarithmic column mass has been chosen 
to correspond to the range in height for y=0. This figure il¬ 
lustrates clearly that the corrugation of the transition region is 
much smaller on a column mass scale than on a height scale 
and that the contribution to intensity, both at line core and for 
the integrated total intensity, is concentrated to a small range 
in column mass. We can also see that indeed the contribution 
to intensity is around a column mass of 10 -6 gcm -2 . 

In sum we expect line core diagnostics (especially the shift 
of the core) to give information from the region just below the 
transition region in the quiet Sun. Diagnostics using the full 


line (like moments of the intensity) will be influenced by a 
much larger part of the chromosphere. We will return to the 
statistical correlation between various line shift diagnostics 
and the physical velocity in the simulation in Paper II. 



0 5 10 15 20 
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Figure 15. Contribution function to intensity at line core (top) and the con¬ 
tribution function to total intensity (bottom) on a height scale along the cut 
through the Bifrost atmosphere at X=12 Mm shown as grey-scale images. 
Temperature contours are given at 40 kK (red) and 10 kK (blue). The in¬ 
tegration for the contribution function to total intensity was performed over 
±20 km s -1 
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Figure 16. As Figure 15 but on a logarithmic column mass scale (see Equa¬ 
tion 6) 


5.2. Line width 

For an optically thin line formed in the CRD regime, the 
width of the intensity profile is directly given by the width of 
the absorption/emission profile 

d> v = -=-—e W (7) 

\Z7lAv D 

where Av is the difference in frequency from the line centre 
frequency, Vo, and 


Av d = ^ 7 — +£2 (8 ) 

c V m 

where ^/2kT/m gives the thermal Doppler broadening and 
£ is a non-thermal velocity. Here k is the Boltzmann con¬ 
stant, m is the mass of the atom/ion and we have assumed a 
Maxwellian velocity distribution at the temperature T. In an 
isothermal plasma and in the absence of a macroscopic veloc¬ 
ity field we get a symmetric Gaussian intensity profile with 
maximum intensity at Vq and a full-width of 2Av D at 1 je of 
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the maximum intensity. Other popular measures of line width 
include the standard deviation, <7, of a Gaussian fit 

_ Av 2 

Iy = {Io-I c )e 2 - 2 +/c, (9) 

the full-width at half the maximum intensity (WfwhmX and 
the second moment of intensity 

w 2 = \j /(v-vi) 2 (/ v -/ c )Jv// tot (10) 

where 

4,t = j(Iv-Ic)dv (11) 

is the zeroth moment and 

V\ — J Av(/ V -I c )dv/I tot (12) 

is the first moment. 

For a Gaussian intensity profile we have the following rela¬ 
tions between the four width measures: 

Av d = V2o (13) 

Wfwhm = 2v / ln2AV/) (14) 

W 2 = o. (15) 

For an optically thick line (as our C II lines) the situation 
is different. The width of the absorption profile still plays a 
role but in addition the source function variation with depth 
is crucial. To illustrate this we have carried out a full 3D 
non-LTE calculation with the Bifrost atmosphere described 
earlier, but with all velocities set to zero. We thus have no ef¬ 
fects from non-thermal velocities and the CII 133.4nm inten¬ 
sity profile is symmetric. In Figure 17 we show the intensity 
profiles and corresponding source functions for two columns 
in this simulation. The red profile is illustrating a narrow, 
single-peaked profile where the source function increases al¬ 
most linearly with decreasing logarithmic column mass from 
the point in the atmosphere where optical depth is unity in 
the continuum (log 10 m c = —2.5) to optical depth unity in the 
core (log 10 m c = —5.6). The temperature in the line-forming 
region is 11 kK which gives a thermal broadening (^/2kT/m) 
of 4 km/s and a thermal Wfwhm = 6.7 km/s. For the red pro¬ 
file we find Wfwhm = 9.3 km/s which is a factor 1.4 larger 
than the thermal width, even in this model where there are 
no non-thermal motions. This is sometimes called opacity 
broadening of an optically thick line. This broadening is de¬ 
pendent on the source function variation with depth. If the 
source function has a peak in the low atmosphere, just above 
the formation height of the continuum, we get steep emission 
flanks, a double peak intensity profile and a larger width. The 
blue profile in Figure 17 is an example of such a scenario with 
Wfwhm = 15.9 km/s with a temperature in the formation re¬ 
gion of 9.6 kK. 

We can get an estimate of the ’’opacity broadening” factor 
by using the Eddington-Barbier relation (that states that the 
emergent vertical intensity is equal to the source function at 
optical depth unity). 

The optical depth is given by 

dz v = -Xvdz (16) 

where Xv is the opacity per volume. The column mass is given 
by 

(17) 




Figure 17. Intensity as function of wavelength (given as equivalent velocity) 
at two columns in the Bifrost atmosphere with zero velocities (left) and cor¬ 
responding source functions as function of logarithmic column mass (right). 
The intensities at the line core are marked with asterisks in the left panel 
and the corresponding height for optical depth unity in the right panel. The 
heights for optical depth unity in the continuum are marked with vertical 
dotted lines. The points corresponding to the FWHM width are marked with 
plus-signs in the left panel with widths of 9.3 km/s (red) and 15.9 km/s (blue). 
The two positions are at (x,y) = (12,7.5)Mm (red) and (x,y) = (12,6.1)Mm 
(blue). 

with p being the mass density. We thus have 

dty = ^dm c (18) 

P 

and we get the column mass where T v = 1 at frequency Av 
from 

rm c (Av) y 

1= / — dm c (19) 

Jo p 

In an isothermal atmosphere with no non-thermal velocities 
we have a constant atomic absorption profile and for a reso¬ 
nance line of the dominant ionization stage we get a constant 
line opacity per unit mass: 

T/v 

^ = aie W (20) 

with ai a constant. We assume the same to be true also for the 
continuum opacity, Xc (where Xv = Xiv+Xc)- We then get 


1 y 

- TT~\ = a l e +- 

m c ( Av) p 


( 21 ) 


Using the same equation for the continuum, Av «», we can 
eliminate the Xc/p term: 


1 


m c ( Av) m c (°°) 

For the line centre, Av = 0, we get 
1 1 


1 -(#M 2 

= aie {Av d j 


m c ( 0) m c (°o) 


= a t 


This gives 


Av 


A Vd \ 


In 


m c (0) m c ( oo) 


1 


1 


m c { Av) m c ( °°) 

In the case m c (°°) m c ( Av) this reduces to 


( 22 ) 


(23) 


(24) 


Av / T m c ( Av) 

a^ = y n ^oT 


dm c = —pdz 


(25) 
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For the blue profile in Figure 17 we have for the FWHM for¬ 
mation height log 10 (m c (Av)) « —2.7, log 10 (m c (0)) = —5.8 
and log 10 (m c (oo)) = —2.4. Using Equation 24, this gives an 
opacity broadening factor of 2.8, which is similar to the actual 
value of 15.9/6.1=2.6. For the red profile we get an opacity 
broadening factor of 1.9, which is larger than the actual value 
of 1.4. This is because the Eddington-Barbier relation under¬ 
estimates the intensity peak when the source function is very 
steep and non-linear in the line forming region and therefore 
overestimates the value of m c (Av) that is necessary to reach 
optical depth unity. Equation 24 also gives the correct po¬ 
sition of the emission peaks of the blue profile. We conclude 
that the width of the CII lines can be 1.4-3 times what is given 
by the width of the atomic absorption profile with larger val¬ 
ues for double peaked profiles. 

5.3. Intensity ratio 

The intensity ratio between the components in the multiplet 
gives additional diagnostic value. Since the two levels within 
each term are very close in energy, collisions will be efficient 
in driving the population ratio towards the ratio of the degen¬ 
eracies. The lower levels have a ratio of their degeneracies 
of two and our simulations confirm that the population densi¬ 
ties have this ratio throughout the line forming regions. The 
opacity scales as nifi u where ni is the population density of 
the lower level and fi u is the absorption oscillator strength. 
The oscillator strengths of the two lines have a ratio of 0.9 
such that the opacity in the 133.5708nm line is a factor of 1.8 
times that of the CII 133.4 nm line. In the optically thin limit, 
we have 

! v- J (Xiv+Xc)S v ds (26) 

In the coronal approximation (collisional excitation, radiative 
deexcitation), the source functions of the two lines are the 
same and the line opacity, Xiv> has the above ratio of 1.8 be¬ 
cause the atomic absorption profiles are also identical. The 
intensity ratio is thus 1.8 at all frequencies across the line 
profiles where the line opacity dominates over the continuum 
opacity. 

For an optically thick line, this is no longer true. We start 
by discussing the case of double peaked profiles and equal 
source functions. The two emission peaks of each line are 
formed at the local source function maximum and we there¬ 
fore get equal peak intensities (but with a larger separation 
for the stronger C II 133.5 nm line). The self reversal is 
formed higher up where the source function is decreasing 
with height and the stronger line has its core formed higher 
and therefore has a lower core intensity. Single peak profiles 
mean that the source function is increasing with height and 
the C II 133.5 nm line will therefore have higher peak inten¬ 
sity than the C II 133.4nm line. In general, the two source 
functions are not equal and we can in principle get any in¬ 
tensity ratio. The weakest blend share upper level with the 
C II 133.4nm line and since it becomes optically thin deeper 
down than the other lines, it tends to depopulate the upper 
level decreasing the source function also of the C II 133.4nm 
line. This source function difference tends to increase the in¬ 
tensity ratio between the lines. 

We illustrate the dependency of the ratio on the source func¬ 
tion behaviour in Figure 18 where we show intensity profiles 
and source functions for the two lines at two locations in the 
Bifrost simulation. The velocities were set to zero before the 
Multi3D calculation of intensities. We thus get symmetric in¬ 


tensity profiles for the C II 133.4nm line but asymmetric pro¬ 
files for the C II 133.5 nm line because of the blend in the blue 
wing. The profiles in red color were taken from the column 
at ( x,y ) = (12,21) Mm where we have single peak intensity 
profiles (left panel). The locations of line core optical depth 
unity are marked on the source functions displayed in the right 
panel. The C II 133.5nm line core (dashed line) is formed 
higher than the C II 133.4nm line (lower column mass) and 
that would by itself give a higher core intensity since the 
source function is increasing with decreasing column mass. 
The main reason for the high intensity ratio of 1.9 is the large 
difference in the source functions. The profiles in blue color 
were taken from the column at (x,y) = (12,6.1) Mm (same as 
the blue colored profile in Figure 17) where the intensity pro¬ 
file has two peaks. For the stronger C II 133.5 nm line we do 
not reach low enough opacity on the blue (short wavelength) 
side of the stronger component to reach optical depth unity 
at the source function maximum depth and we need to go all 
the way to the blue side of the blend to get the corresponding 
intensity peak. This explains the rather weird profile (dashed 
blue line in the left panel). The red (long wavelength) peaks 
have an intensity ratio of 1.4, again mainly caused by differ¬ 
ences in the source functions. 

The source function differences are driven by radiative 
rates, both losses in the weakest line decreasing the source 
function of the C II 133.4nm line that shares the upper level 
with the weak blend and absorptions of these photons in 
the C II 133.5 nm line (which lead to an increase of the 
C II 133.5 nm source function). Collisions between the upper 
2 p 2 D levels are not efficient enough to drive the source func¬ 
tions to equality. We have also included collisions with neu¬ 
tral hydrogen between these levels but a factor of 10 higher 
total collisional rates would be needed to significantly affect 
the source function ratios. 

We conclude that the intensity ratio of the C II 133.5 nm 
and C II 133.4nm lines is primarily indicating the ratio of the 
source functions of the lines. This ratio is normally larger than 
one but can have any value, including the optically thin value 
of 1.8. We expect a lower ratio when densities are high in the 
line forming region, as is the case for the intensity peaks of 
double peak profiles. 




Figure 18. Intensity as function of wavelength (given as equivalent veloc¬ 
ity) at two columns in the Bifrost atmosphere with zero velocities (left) for 
the C II 133.4nm line (solid) and the C II 133.5nm line (dashed) and cor¬ 
responding source functions as function of logarithmic column mass (right). 
Maximum intensities are marked with asterisks in the left panel and the cor¬ 
responding height for optical depth unity in the right panel. The two Bifrost 
columns are at (x,y) = (12,21)Mm (red) and (jc,y) = (12,6.1)Mm (blue). 
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5.4. Comparison with optically thin modelling 

We have found that the C II lines are mostly formed in the 
optically thick regime. If we instead assume optically thin 
formation with the ionization equilibrium set by collisional 
ionization and dielectric recombination without density de¬ 
pendent effects (standard coronal approximation) we severely 
underestimate the amount of C II in the upper chromosphere 
and overestimate the amount in the higher temperature region 
(see Figure 7). In the optically thin regime, the intensity is set 
by the source function multiplied by the optical thickness of 
C II in the region. With the low densities where C II exists in 
the coronal approximation, this gives an optical thickness of 
less than 0.1. The source function then needs to be more than 
ten times higher than what it is in our model at optical depth 
unity to give a similar intensity. This is normally not the case 
and we thus get a weaker line from the coronal approximation 
than in our non-LTE modelling. This is shown in Figure 19 
where we compare the emergent intensity from assuming the 
coronal approximation (left panel) with our simulation em¬ 
ploying the full 3D non-LTE modelling (right panel). 



x[Mm] x[Mm] 

Figure 19. Emergent intensity at 133.4663 nm from the coronal approxima¬ 
tion (left panel) and from the full 3D non-LTE calculation (right panel). 

6. SUNSPOT MODEL 

In addition to the quiet sun case, we have also studied the 
behavior of the lines in the sunspot model atmosphere of Lites 
& Skumanich (1982). The temperature structure of this model 
is shown in Figure 2. The photospheric temperature is natu¬ 
rally lower than in the quiet Sun model VAL3C and a low tem¬ 
perature of 3.5 kK continues up to logj 0 m c = — 2 where there 
is a rapid temperature increase to a chromospheric plateau 
of 6.5 kK, slightly above the VAL3C temperatures at corre¬ 
sponding column masses. There is a rapid temperature in¬ 
crease in the transition region at log 10 m c = —5.44. A temper¬ 
ature plateau has been introduced at a temperature of 24 kK 
in order to get a Lyman-a flux similar to the observed flux. 

The emergent intensity for the C II 133.5 nm lines is shown 
in Figure 20. The ratio of the maximum intensities of the 
two lines is 1.6. The line profiles deviate from a Gaussian 
shape in that the core part is flatter (especially for the stronger 
C II 133.5 nm line) and the flanks are steeper than the best- 
fit Gaussian. The FWHM of the intensity profile of the un¬ 
blended CII 133.4nm line is 16.4 km s _1 . The temperature at 
the formation height is 24 kK and the microturbulence in the 
sunspot model is 2.6 km s _1 . In an optically thin formation, 
the thermal width and the microturbulence would combine to 
give a FWHM of 10.5 km s _1 . We thus have an opacity broad¬ 
ening factor of 1.6. 

The intensity formation of the C II 133.4nm line in this at¬ 
mosphere is shown in Figure 21. The velocities for this atmo- 



Figure 20. Intensity as function of wavelength (in equivalent velocity units) 
for the C II 133.4nm line (solid) and the C II 133.5nm line (dashed) in the 
sunspot umbral model of Lites & Skumanich (1982). 


sphere are set to zero. The T v = 1 height (red dashed line in 
all panels) ranges from 0.8 Mm in the continuum to 2.0 Mm 
at the line center. The source function (solid green in upper 
right panel) is equal to the Planck function (yellow dashed line 
in the upper right panel) from the photosphere up to 0.9 Mm 
height where it decouples. The source function continuously 
increase with the increase in temperature and becomes max¬ 
imum in the transition region. The monotonically increasing 
source function leads to a single peak emission line (white 
curve in the lower right panel). The contribution function to 
intensity (image in lower right panel) is concentrated around 
2.1 Mm height around optical depth unity which is also the 
lower part of the plateau at 24 kK in the model. Much of the 
line intensity is thus formed at the temperature plateau which 
accounts for the rather flat line core. This is even more pro¬ 
nounced for the stronger CII 133.5 nm line where a significant 
fraction of the line core is formed at the temperature plateau 
where the source function is rather constant (upper left panel). 

7. DISCUSSION AND CONCLUSIONS 

The ionization balance and intensity profiles of the C II lines 
at 133.4-133.6 nm can be modelled with a 9-level model atom 
assuming complete frequency redistribution (CRD). The ion¬ 
ization balance C I/C II is dominated by photoionization and 
radiative recombination. The balance C II/C III is dominated 
by collisional ionization and dielectronic recombination. In 
the VAL3C atmosphere, carbon is more than 50% in the form 
of C II in the temperature range 5.6-35 kK, where especially 
the low temperature end is atmosphere dependent since it is 
set by photoionization which is a non-local process. Hori¬ 
zontal scattering is important for the line core intensity which 
necessitates 3D non-LTE modelling and makes line core di¬ 
agnostics influenced by non-local conditions. 

The line core is formed at temperatures 6-25 kK at a col¬ 
umn mass close to 10 -6 g cm -2 . The continuum formation is 
at 300-800 km and is dominated by background opacity from 
neutral silicon. This background opacity needs to be treated 
in non-LTE in order to get the continuum intensity right but 
assuming LTE only affects the continuum intensity and not 
the line core intensity. The lower continuum formation height 
is for kilo-Gauss flux concentrations in the photosphere. 

The atomic absorption profile width is set by thermal and 
non-thermal broadening. The emergent intensity profile is 
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Figure 21. As Figure 11 but for the sunspot model atmosphere of Lites & 
Skumanich (1982). 

broader than the atomic absorption profile because of the op¬ 
tical thick formation. This excess ’’opacity broadening” is 
a factor of 1.2-4. The smaller opacity broadening happens 
for single peak intensity profiles that originate from a source 
function with a steep increase into the transition region and 
with a low chromospheric temperature. The larger opacity 
broadening is for double peak intensity profiles with central 
reversals that result from a source function with a local max¬ 
imum in the lower chromosphere because of a steep tempera¬ 
ture rise there. 

The intensity ratio between the C II 133.5 nm and the 
C II 133.4nm line is lower than the optically thin value of 
1.8 in the VAL3C model, the umbral model as well as for 
most columns in the Bifrost snapshot. The optically thick 
formation can give any ratio, depending on the ratio of the 
source functions of the two lines. The ratio is lower for dou¬ 
ble peak profiles where the intensity peaks are formed in the 
lower chromosphere than for single peak profiles formed in 
the transition region. 

The C II lines are formed in the optically thick regime and 
erroneously assuming an optically thin formation using the 
coronal approximation (neglecting photoionization) leads to 
too high formation temperatures and too low intensities. 

In the next paper in this series we will discuss statistical 
correlations between atmospheric parameters and observables 
and explore the diagnostic potential of the C II lines. 
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